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Abstract 

The Mars Global Surveyor (MGS) arrived at Mars in September 1997 near Mars’s southern 
spring equinox and has now provided monitoring of conditions in the Mars atmosphere for more 
than half a Mars year. The large majority of the spectra taken by the Thermal Emission 
Spectrometer (TES) are in a nadir geometry (downward looking mode) where Mars is observed 
through the atmosphere. Most of these contain the distinct spectral signature of atmospheric 
dust. For these nadir-geometry spectra we retrieve column-integrated infrared aerosol (dust) 
opacities. TES observations during the aerobraking and science-phasing portions of the MGS 
mission cover the seasonal range L s = 184°-28°. Excellent spatial coverage was obtained in the 
southern hemisphere. Northern hemisphere coverage is generally limited to narrow strips taken 
during the periapsis pass but is still very valuable. At the beginning of the mission the 9-/rm 
dust opacity at midsouthern latitudes was low (0.15-0.25). As the season advanced through 
southern spring and into summer, TES observed several regional dust storms (including the 
Noachis dust storm of November 1997) where peak 9-/im dust opacities approached or exceeded 
unity, as well as numerous smaller local storms. Both large and small dust storms exhibited 
significant changes in both spatial coverage and intensity over a timescale of a day. Throughout 
southern spring and summer the region at the edge of the retreating southern seasonal polar ice 
cap was observed to be consistently more dusty than other latitudes. 
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1. Introduction 

Dust aerosols are always present in the Mars atmo- 
sphere. Their presence significantly effects the ther- 
mal structure of the atmosphere and is a major driver 
of atmospheric circulations at all spatial scales. A 
long history of both Earth-based and spacecraft-based 
observation exists of dust in the Mars atmosphere. 
Earth-based observations of Martian dust date back 
at least to the 1870s with the work of Schiaparelli. 
Martin and Zurek [1993] give a good review of the his- 
tory of Earth-based (and spacecraft-based) dust cloud 
observations. More recently, images from the Hub- 
ble Space Telescope [e.g., James et at. , 1997; Wolff et 
al 1997] and observations from amateur astronomers 
[Beish and Parker , 1990] have helped to character- 
ize the frequency and location of dust storms. Al- 
though clouds of dust can be easily identified in im- 
ages of Mars, it is very difficult to accurately estimate 
a column-integrated opacity. 

Spacecraft observation of Mars has allowed a more 
quantitative analysis of dust opacity. Data from the 
infrared spectrometer carried by Mariner 6 and 7 
clearly show the signature of dust at 900-1200 cm -1 
(8-11 //m) [Horn et al., 1972: Cimino and Calvin , 
1996]. Combined with an estimate for atmospheric 
and surface temperatures, these data can be used to 
compute column-integrated dust opacity using radia- 
tive transfer techniques. Mariner 9 also carried an 
infrared spectrometer [Hanel et al .. 1972a] and mon- 
itored the decay of the global dust storm of 1971 
[Hanel et al 1972b; Toon and Pollack , 1977]. Dust 
opacities have recently been derived from these data 
by Santee and Crisp [1993] and by Fenton et al. 
[1997]. The two Viking orbiters did not carry infrared 
spectrometers, but data from the broadband infrared 
thermal mapper (IRTM) [Chase et al., 1978] can be 
used to accurately retrieve dust opacity and have 
proven to be very useful for monitoring dust opaci- 
ties and local dust storms during the Viking mission 
[ Martin , 1986; Martin and Richardson , 1993; Peter - 
freund , 1985]. Observations of the Sun by the Viking 
lander give dust opacity at visible wavelengths [Pol- 
lack et al., 1979; Colburn et al.. 1989]. Thermal in- 
frared spectra recorded using ground-based telescopes 
[Pollack et al ., 1990] and the Kuiper Airborne Obser- 
vatory [ Roush et a/., 1995] have also been used to es- 
timate dust opacities. Most recently, the Pathfinder 
lander has used observations of the Sun during the 
day and of bright stars at night to obtain dust opac- 
ity [Smith et al., 1997; Smith and Lemmon , 1999]. 


In this paper we extend this baseline of observa- 
tions by reposing results of the retrieval of column- 
integrated dust opacity using infrared spectra taken 
by the Thermal Emission Spectrometer (TES) on- 
board the Mars Global Surveyor (MGS). At least two 
other instruments on MGS can also provide useful 
information on dust opacities. The Mars Observer 
Camera (MOC) imaged the large regional dust storm 
in Noachis Terra, as well as numerous other smaller 
storms [e.g., Malm et al., 1998a; Matin, 1998], and 
Toigo et al. [1998] have used an analysis of shad- 
ows in MOC images to estimate visible dust opacities. 
Ivanov and Muhleman [1998] have used data from the 
Mars Orbiter Laser Altimeter (MOLA) to infer dust 
opacity from measurements of atmospheric transmit- 
tance. 

We next review the characteristics of the TES in- 
strument and its observations in section 2. In sec- 
tion 3 we describe the retrieval algorithm used to ob- 
tain column-integrated dust opacities. In section 4 we 
present results and compare them with results from 
previous spacecraft missions. Finally, we summarize 
our findings in section 5. 


2. Data Set 


2.1. TES Instrument 


The Thermal Emission Spectrometer (TES) is a 
thermal infrared interferometer/ spectrometer with 
additional broadband visible and thermal channels 
[Christensen et al., 1992]. Six detectors in a three- 
by-two array simultaneously take spectra covering the 
spectral range from 200 to 1600 cm -1 (6-50 ^m), 
with a selectable sampling of either 5 or 10 cm -1 (a 
large majority of the data collected so far has 10- 
cm -1 sampling). A pointing mirror allows TES to 
view from nadir to above both the forward and aft 
limbs, where the atmosphere is observed without di- 
rect contribution from the surface. Each pixel sub- 
tends a 8.3-mrad field of view. The highly eccentric 
orbit of MGS during the period of observations con- 
sidered here causes highly variable spatial resolution 
for TES observations, ranging from 10 km (including 
smear caused by spacecraft motion) near periapsis to 
several hundred kilometers near apoapsis. 

Figure 1 shows examples of typical daytime TES 
spectra. Spectra are shown here in terms of their 
equivalent brightness temperature to emphasize fea- 
tures. The C0 2 gas absorption centered at 667 cm" 1 
(15 pm) is used to retrieve atmospheric temperatures 
[see Conrath et al., this issue]. The spectral signa- 
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tures of dust and H 2 O ice aerosols are clearly visible 
and are distinct from each other and from the 15- 
/im CO 2 band. Dust has a broad peak centered at 
1075 cm" 1 and nonnegligible absorption throughout 
the entire spectrum except near 1300 cm -1 . Absorp- 
tion from HoO ice has a broad peak near 825 cm" 1 
and a sharper peak at 229 cm" 1 (only hinted at in 
Figure 1) and is small for wavenumbers > 1000 cm" 1 
[see Curran et a/., 1973; J. C. Pearl, Mars water ice 
clouds: Observations by the Thermal Emission Spec- 
trometer (TES) during the first Martian year, submit- 
ted to Journal of Geophysical Research, 1999]. The 
most transparent part of the spectrum is near 1300 
cm -1 . 

2.2. Mars Global Surveyor Orbit 
Characteristics 

The Mars Global Surveyor (MGS) entered into or- 
bit around Mars in September 1997 and performed 
over 200 orbits of aerobraking before entering a park- 
ing, or “science-phasing,” orbit in March 1998. Mars 
Global Surveyor remained in the science-phasing or- 
bit until September 1998, when aerobraking to the 
final mapping orbit resumed. A brief summary of the 
first part of the mission is given by Albee et al. [1998]. 

The initial orbital period of MGS was ^~45 hours. 
During the first period of aerobraking (September 
1997 to March 1998) the orbital period was gradu- 
ally decreased to 11.6 hours, the latitude of periapsis 
moved from 32°N to 51°N, and the local time of peri- 
apsis changed from 1820 hours to 1220 hours. There 
was no aerobraking attempted during the science- 
phasing orbit, and so the orbital period remained at 
11.6 hours. The latitude of periapsis continued mi- 
grating northward from 51°N going over the north 
pole and reaching 61°N at the end of the science- 
phasing observations. Local time of periapsis was 
near dawn at the beginning of science phasing and 
then near dusk after periapsis moved over the north 
pole. 

2.3. Observations During Aerobraking and 
Science Phasing 

Because of the elliptical nature of the MGS orbit, 
observations from TES fall into two main categories: 
those taken during the periapsis pass and those taken 
during the rest of the orbit (the “apoapsis portion” 
of the orbit). Observations were split fairly evenly 
between the two categories during aerobraking when 
the orbital period was long. Data playback was lim- 
ited during the science-phasing portion of the mission, 


and the majority of observations were focused on the 
periapsis pass. 

For any one orbit the periapsis pass gives one nar- 
row (three-pixel-w’ide) strip of observations running 
roughly north to south. Coverage was primarily in 
the northern hemisphere and was usually in daylight. 
During the apoapsis portion of the orbit the space- 
craft was placed into a roll configuration so that by 
using the TES pointing mirror, most or all of the disk 
of Mars could be imaged once per 100 min. This al- 
lowed large portions of the southern hemisphere to be 
sampled, although more than half of the spectra were 
taken on the nighttime portion of the planet. 

In this paper we use spectra with surface tem- 
peratures warm enough (> 220 K) to provide suffi- 
cient thermal contrast between the surface and at- 
mosphere for the accurate retrieval of atmospheric 
opacities. This restricts coverage to daytime measure- 
ments. Spectra from both the periapsis pass (which 
was in daylight during aerobraking and science phas- 
ing) and the daytime portions of the apoapsis observa- 
tions are used. When daytime periapsis and apoapsis 
data overlap spatially and are separated in time by 
less than a day, dust opacities between the two usu- 
ally agree to within measurement error. 

3. Dust Opacity Retrieval Algorithm 

To retrieve dust opacity, we first compute the 
equivalent column-integrated opacity of pure absorbers 
as a function of wavenumber. Then, we estimate the 
contribution of dust to the total opacity by fitting 
predetermined spectral shapes (opacity as a function 
wavenumber) for dust, water ice, and the effect of 
a nonunit emissivity surface to the observed opacity 
spectrum. 

3.1. Opacity Spectrum 

If we neglect the small contributions from the so- 
lar beam and scattering from atmospheric aerosols, 
and assume a plane-parallel atmosphere, then the 
observed monochromatic radiance of Mars a s a fre- 
quency of wavenumber, / 0 bs(^), can be written as 

W") = <(") B[T sur{ ,u]e- T ^ fi 



where e(u) is the surface emissivity at frequency v, r 0 
is the normal column-integrated aerosol opacity, p is 
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the cosine of the emission angle, B[T . v] is the Planck 
function, T sur f is the surface temperature, T(r) is 
the atmospheric temperature, and the integral is per- 
formed from the spacecraft (at r = 0) to the surface 
(at r = r 0 ). 

The procedure used to retrieve column-integrated 
aerosol opacity is a numerical solution of (1) for tq(u). 
Input are the data (radiance at each wavenumber), 
a priori information in the form of a previously re- 
trieved surface temperature and atmospheric temper- 
ature profile, and assumptions that the dust is well- 
mixed and nonscattering that are described in detail 
below. The expected radiance for a range of opacities 
(r 0 ) is computed using (I) until a pair of opacities is 
found that brackets the observed radiance. The opac- 
ity where the expected radiance matches the observed 
radiance is then found by standard root-finding tech- 
niques. This procedure is repeated for each wavenum- 
ber in the observed spectrum. The final solution is 
the spectral dependence of the equivalent column- 
integrated opacity of purely absorbing aerosols for a 
ray normal to the surface, ro(t')- 

To solve (1) for opacity, atmospheric and surface 
temperatures are required. Atmospheric tempera- 
tures were retrieved for each TES spectrum using a 
constrained linear inversion of radiances in the 15-/im 
CO 2 band (see the accompanying paper by Conrath 
et ai . [this issue] for details). For this study the re- 
trievals for temperature and aerosol opacity were per- 
formed sequentially (instead of simultaneously). In 
the temperature retrieval we have accounted for dust 
opacity to first order by using an effective surface tem- 
perature calculated by taking the brightness tempera- 
tures in narrow spectral intervals on either side of the 
COo band and averaging them. The spectral intervals 
are 508-529 and 814-825 cm” 1 . Temperatures from 
the surface to about 0.1 mbar (35 km) are obtained 
with a typical vertical resolution of about one scale 
height (10 km). Surface temperatures were obtained 
by using the brightness temperature in a small spec- 
tral window near 1300 cm” 1 , where the brightness 
temperature is highest (see Figure 1). 

We make two major simplifying assumptions to 
be able to numerically solve (1) for the equivalent 
column-integrated opacity tq. First, we assume that 
atmospheric opacity sources (primarily dust and wa- 
ter ice) are well-mixed with the CO 2 gas. Obviously, 
this is not always the case. Using Viking orbiter 
images, Jaquin et ai [1986] observed vertical struc- 
ture in dust hazes. However, analysis of limb scans 
from the Viking lander Sun diodes [Pollack et ai , 


1977] and images from the Pathfinder lander [Smith 
et a/., 1997] indicates that at most times the well- 
mixed approximation is a good one. Furthermore, 
our experience with the TES data shows that dust 
opacity scales very well with both surface pressure 
and l/cos(emission angle), both indicating that dust 
opacity is proportional to the amount of gas along 
the line of sight. Because well-mixed dust means that 
dust opacity along the line of sight is directly propor- 
tional to atmospheric pressure, the integration over r 
in (1) can be changed to an integration over pressure. 
Atmospheric temperature T(r) can then be replaced 
with T(p), which is the quantity that is directly re- 
trieved from the 15-/im CO 2 band. 

The second assumption that we make is that the 
dust aerosols are nonscattering. Therefore we are re- 
trieving an equivalent absorption dust opacity instead 
of the full extinction (absorption plus scattering) dust 
opacity. For nadir viewing geometry and the typical 
opacities found outside of dust storms of r 0 < 0.25, 
the effect of scattering is small, and the absorption 
dust opacity is very nearly equal to the extinction 
dust opacity. For the daytime spectra used here, the 
absorption dust opacity tends to be a slight underes- 
timate (typically 10-20%) of the true extinction dust 
opacity. During the very dustiest times when opac- 
ities approach unity, scattering may become impor- 
tant, and in that case our dust opacity estimates be- 
come more of an indicator of relative dust opacity 
rather than an absolute measure of the amount of 
dust in the atmosphere. 

3.2. Spectral Shape Fitting 

In order to obtain an accurate dust opacity, it is im- 
portant to account for the opacity contribution from 
water ice and from a nonunit emissivity surface, espe- 
cially during times of relatively low dust opacity. Not 
including these contributions can lead to a substantial 
overestimation of dust opacity. A detailed descrip- 
tion of this effect is given by Smith et ai [this issue]. 
To estimate the contributions of dust, water ice, and 
nonunit surface emissivity to the opacity spectrum 
found above, we assume that the opacity spectrum 
t 0 (i/) can be described by 

tq(is) — - A, j ust (lat , lon)/dust (^) Ai ce (lat., lon)/i ce (i/) 

+ A sur f (lat, lon)^ sur f(^) — , (2) 

where /d ust and f\ ce are the spectral shapes (opac- 
ity as a function of wavenumber) of dust and water 
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ice aerosols, < 7 sur f is the spectral shape of the surface 
(one minus emissivity as a function of wavenumber), 
the .4 parameters describe the horizontal (spectrum- 
to-spectrum) variation in the contribution from each 
component, and t is surface emissivity. 

For each opacity spectrum a least squares fit for 
Adust, Ace* and -4 surf is performed by assuming that 
the spectral shapes /dust, /ice, and g suv f are fixed and 
known [see Smith et a/., this issue; Bandfield et a/., 
this issue]. The derivative term in (2) is computed nu- 
merically for each spectrum by varying surface emis- 
sivity in (1). Finally, a simple quality control is per- 
formed that excludes all fits with high residual and 
unphysical values for the A parameters. The resulting 
best fit value for Adust is then the column-integrated 
opacity of dust at 1075 cm -1 (since /dust = 1 at 1075 
cm -1 ) and is the quantity for "dust opacity” that we 
use in the remainder of this paper. 

In performing the least squares fit, we use the spec- 
tral shapes for dust and water ice opacity adopted by 
Smith et al. [this issue]. For the surface emissivity 
spectral shape we use that derived for orbit 219 (Cim- 
meria Terra) by Smith et al. [this issue], which has a 
broad feature from 850 to 1200 cm" 1 and an emissiv- 
ity at 250-550 cm -1 roughly equal to that at 1000- 
1100 cm” 1 . Although the surface spectral shape has 
some similarity to the dust spectral shape at 850-1200 
cm” 1 , the opacity of dust is much reduced at 250— 
550 cm” 1 relative to the peak at 1075 cm” 1 , while 
the surface emissivity has comparable values in the 
two spectral regions. Surface emissivity varies with 
location on the planet but is always > 0.90 at 1000 
cm -1 . 

3.3. Estimation of Uncertainties 

Uncertainty in the retrieved value of dust opac- 
ity comes from instrument noise and calibration, er- 
rors in the surface temperature and atmospheric tem- 
perature, and the assumptions outlined above (fixed 
spectral shapes and a well-mixed vertical distribution 
for dust). Using a conservative (noisy) estimate for 
the instrument noise of 3 x 10” 8 W cm” 2 sr” l /cm” 1 
[Christensen et al., 1992, 1998: Christensen , 1999] 
gives a typical dust opacity uncertainty of 0.02 for 
one spectral channel. However, because we fit spectral 
shapes over many frequencies, the uncertainty caused 
by noise is very small. Conservative estimates for the 
uncertainty in surface temperature (0.5 K) and in at- 
mospheric temperatures (3 K) also give dust opacity 
uncertainties of 0.02 each. All of these uncertainties 
are inversely related to the amount of thermal con- 


trast between the atmosphere and surface. Values 
quoted here are for a thermal contrast of 20-40 K, 
which is. again, conservative for the observations used 
in this study. 

The uncertainty introduced by assuming fixed spec- 
tral shapes is much more difficult to quantify. Band- 
field et al. [this issue] and Smith et al. [this issue] 
show that dust and water ice spectral shapes have 
only minor variations (< 5%) for the vast majority of 
spectra. The surface emissivity spectral shape may 
show more variation. There are certainly regions that 
display surface emissivities that are significantly dif- 
ferent than the spectral shape computed over Cimme- 
ria Terra on orbit 219 (for example, the region with 
strong hematite features described by Christensen et 
al. [this issue]. However, Christensen [1982, 1998] 
shows that the presence of surface emissivity spectral 
features is highly correlated with low-albedo regions 
and that although there are meaningful differences 
between surface emissivity spectra taken at differ- 
ent low-albedo locations, the general spectral shape 
of those surface emissivity spectra is broadly similar 
to the orbit 219 Cimmeria Terra spectrum. Further- 
more, any surface emissivity that is significantly dif- 
ferent than what we assumed here would cause a high 
residual in our fit and would be rejected by our quality 
control procedure. Thus, for the purpose of obtaining 
dust opacity, most surface emissivity spectral shapes 
computed thus far can be considered to be variations 
of the orbit 219 shape. We estimate that the uncer- 
tainty in dust opacity introduced by assuming a fixed 
surface emissivity spectral shape is no more than 0.03. 

The uncertainty caused by dust not being well- 
mixed (dust layering) is also very difficult to estimate. 
However, as noted earlier, our experience with the 
data is that this effect is small and has negligible im- 
pact on the uncertainty calculation. Adding all these 
uncertainties together, assuming that they are uncor- 
related, leads to a total uncertainty estimate of about 
0.05 for any one dust opacity retrieval. 

4. Results 

The procedure described above yields ^330,000 in- 
dividual retrievals of column-integrated dust opacity. 
The dust opacities presented below are given for 1075 
cm” 1 , which is at the peak of the dust absorption. 
The spectral dependence of dust is shown in Figure 
2 . 

Because the approximation that dust is well-mixed 
with the COo gas is quite good, the large topographic 
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variations on Mars make comparisons of dust opacity 
between one place and another difficult unless the to- 
pography is ^removed*’ by scaling to a reference pres- 
sure. In the plates below we scale all dust opacities 
to an equivalent 6.1-mbar pressure. This scaling also 
corrects for the bias in opacities caused by the annual 
variation of surface pressure, which has an amplitude 
of about 20% [Hess et al ., 1980]. 

4.1. Overview 

Plate 1 gives an overview of dust opacity as ob- 
served by TES. Dust opacities at 1075 cm” 1 scaled 
to an equivalent 6.1-mbar surface are given for seven 
different latitude bands. Dust opacity averages are 
obtained by first placing all the data from one orbit 
into bins 15° wide in latitude and 30° wide in longi- 
tude. Then, the dust opacity for each orbit and lati- 
tude band is obtained by averaging all longitude bins 
for that orbit and latitude that contain data. This 
two-step average is performed to minimize bias caused 
by uneven longitude coverage because TES takes ap- 
proximately the same number of measurements near 
the longitude of periapsis as at all other longitudes 
combined. 

When TES first began observations just after south- 
ern spring equinox, there was a general latitudinal 
gradient in dust opacity with higher opacities to the 
south. A moderate regional dust storm at mid- 
southern latitudes was observed at around L a = 190°. 
Throughout the period L, = 184 o -220° there was in- 
creased dust activity at the edge of the retreating 
southern seasonal polar ice cap (see kieffer et ai. [this 
issue] for a detailed description of the recession of the 
south polar ice cap). During the period L*=215°- 
225° there was a general trend of increasing dustiness 
in the southern hemisphere perhaps caused by the 
numerous local storms observed along the edge of the 
south polar ice cap. Leovy et al. [1973] postulated 
that local dust storms along the retreating southern 
polar ice cap would contribute to a background opac- 
ity in the southern hemisphere that could make condi- 
tions favorable for the initiation of large dust storms. 

At L a =225°, TES observed a very large regional 
dust storm which originated in the Noachis Terra 
region. Dust opacity increased at all latitudes and 
rapidly spread east and west over the course of the 
next week. However, it is clear from Plate 1 that the 
increase in dust opacity was far greater to the south of 
30° S than in was to the north and that this large re- 
gional storm did not attain truly global status. The 
main body of the dust storm began to subside af- 


ter a few sols, although significant activity remained 
throughout the rest of the aerobraking phase at high 
southerly latitudes (60°-75°S). 

When the science-phasing orbit began at £,=304°, 
the atmosphere was relatively clear. Shortly after 
that, at Lj=309°, TES observed a moderately large 
regional dust storm that originated north of the Ar- 
gyre basin and spread to the south and east. The 
last part of the science-phasing orbit (L J =0°-28°) was 
characterized by gradual clearing from a regional dust 
storm observed in microwave temperature soundings 
by Clancy et ai [this issue] at L,= 340°, punctuated 
by occasional small localized dust storms. 

The variation of atmospheric temperature with 
time over this period (presented by Conrath et al. 
[this issue], Figure 2) show's correlations with the dust 
opacity trend. There are rapid increases of 10-15 K 
in atmospheric temperatures at the 0.3-mbar pressure 
level (about 30-km altitude) in both the northern and 
southern hemispheres associated with the larger dust 
storms mentioned above (the Noachis dust storm at 
L,= 225° and the regional dust storm at L s - 309°). 
The rapid and large amplitude response of the atmo- 
sphere in both hemispheres implies a strong interac- 
tion between dust opacity and the dynamical state of 
the atmosphere. 

4.2. Noachis Dust Storm 

The major dust event during these observations 
was the Noachis dust storm which began on Novem- 
ber 26, 1997 [L s ~ 225°). Plate 2 shows a map incor- 
porating each dust opacity measurement taken during 
the main part of the storm (T a — 225° -233°). The 12 
strips of data running roughly north and south from 
30° N to 30°S represent the periapsis passes from the 
12 orbits used. The even distribution of observations 
at all longitudes southward of 30°S latitude comes 
from the apoapsis portion of the orbit. Our selection 
of data with a warm surface and a significant thermal 
contrast between the surface and atmosphere excludes 
observations taken over the polar ice cap, which at 
this time extended to about 68°S. 

Plate 3 shows the evolution of the storm. The or- 
bital period at this time was about 30 hours. The 
main core of the storm was located in Noachis Terra 
between 20°-55°S latitude and 330°-10°W longitude. 
During the storm, high dust opacities (> 0.5) spread 
to all longitudes poleward of 40°S. The storm prefer- 
entially spread to the south and east. The evolution 
of the storm is qualitatively similar to the regional 
dust storm modeled by Murphy et al. [1995, Figure 
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3] with expansion of the dust first poleward in direc- 
tion and then quickly spreading to all longitudes at 
high southerly latitudes. 

In the northern hemisphere, dust opacities were 
not significantly elevated above prestorm levels until 
orbit 53 (£,= 226°, approximately 4 sols after the be- 
ginning of the storm). Opacities remained elevated 
until at least orbit 57 (L a = 231°. 9 sols after the be- 
ginning of the storm) before a gradual clearing began. 
In the north, there was somewhat more dust opacity 
at longitudes near the storm center (330°-10°W) than 
in the opposite hemisphere (150°- 190° YV), but there 
was considerably less longitudinal structure than in 
the southern hemisphere. The sparse longitude cov- 
erage in the northern hemisphere for each orbit (data 
are primarily from periapsis passes only) limits de- 
tailed analysis of northern hemisphere data. 

In the southern hemisphere, large changes in the 
intensity and areal extent of the dust storm were ob- 
served from orbit to orbit. The main core of the 
Noachis storm quickly intensified, reaching a maxi- 
mum strength at about orbit 53 (£,= 226°, 4 sols af- 
ter the beginning of the storm). At this time the 
main core covered the area from about 20° S lati- 
tude to at least the edge of the polar cap (68°S) and 
from 270°W longitude through the prime meridian to 
60°W longitude. Large gradients in opacity existed 
in the growth phase of the dust storm with opacities 
at and exceeding unity inside the dust storm but still 
at nearly prestorm levels (0.3-0. 4) just outside the 
dust storm. At the height of the Noachis dust storm, 
atmospheric temperatures inside the dust storm were 
10-20 K warmer at the 0.3-mbar pressure level (about 
30-km altitude) than they were at the same latitude 
outside the dust storm. 

Over the next 4 to 5 sols the main core of the 
Noachis core began to slowly dissipate. However, dust 
continued to spread, and opacities of at least 0.5 cov- 
ered all longitudes poleward of 40°S latitude by or- 
bit 57 (7, =231°, 9 sols after the beginning of the 
storm). By orbits 59 and 60 (£,=233°, 12 sols after 
the beginning of the storm), the Noachis storm cen- 
ter was no longer well organized, and a new center 
of activity at 60°-70°S and 180°-270°W had estab- 
lished itself. The Noachis storm center dissipated by 
orbit 70 (£,=241°, 25 sols after the beginning of the 
storm), but the satellite storm at 60°-70°S persisted 
for quite some time, maintaining a higher dust opac- 
ity than surrounding longitudes until about £,=270° 
(see Plate 1). There was significantly more dust opac- 
ity at 180°-270° W longitude than at other longitudes. 


Although the Noachis storm center was gone by or- 
bit 70, dust opacity levels throughout the southern 
hemisphere remained significantly higher than before 
the storm until the gradual clearing of the atmosphere 
brought them back down to prestorm levels at around 
£,=270°. 

By fitting the dust opacity trends in Plate 1 for 
the latitude bands 45°-60°S and 60°-75°S to an ex- 
ponential function, we estimate an exponential decay 
constant of 49 sols for the main Noachis storm. The 
30°-45°S latitude band has a decay constant of 75 
sols, and the decay constant for more northerly lati- 
tude bands (far outside the core of the storm) is even 
longer. 

4.3. Other Dust Storms 

In addition to the Noachis and its satellite dust 
storm, TES observed several other significant dust 
storms. At £ a =309° a compact but intense dust 
storm started at 20°-50°S latitude, 20°-120°W longi- 
tude (the area north and northwest of Argyre). Plate 
4 is a map showing all dust opacity measurements 
obtained during the storm (£ 5 =309°-319°). As the 
storm evolved it spread south and east forming a new 
core of activity at 65°-80°S latitude, 180°-270° W lon- 
gitude that outlived the activity north of Argyre. In 
this way the evolution of the dust storm was remark- 
ably similar to that of the Noachis dust storm. How- 
ever, unlike the Noachis storm the £j~309° storm had 
very little influence outside its core regions. As shown 
in Plate 4, dust opacity remained at low prestorm lev- 
els (0.15-0.20) at all latitudes north of 65°S between 
150°-320°W longitude. Even in the longitude where 
the storm was most active (0°-120°W), there was no 
significant change in northern hemisphere dust opaci- 
ties outside of a couple isolated local dust storms that 
were probably unrelated to the southern hemisphere 
regional storm. 

A smaller and less intense regional dust storm was 
observed by TES at the beginning of our observations 
at £, = 184°. Plate 5 is a map showing all observa- 
tions between L s ~184°-195°. There are two centers 
of increased dust opacity at southern midlatitudes. 
The main one was located at 40°-60°S latitude, 300°- 
340° W longitude (west of Hellas), and a smaller one 
was located at 40°-60°S latitude, 20°-50°VV longitude 
(Argyre). The edge of the south polar ice cap was at 
about 60 S at this time. These two regions were ac- 
tive simultaneously, and both peaked in intensity near 
£, = 192°. Neither of these storms migrated. 

Numerous smaller dust storms were also recorded 
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by TES, ranging in size from several hundreds of kilo- 
meters down to a single TES pixel (roughly 10 km). 
Many of these storms occurred at the edge of the re- 
treating south polar seasonal ice cap. Few of these 
local storms had normalized dust opacities > 0.5, and 
none lasted for longer than a sol or two. Table 1 lists 
a few of the larger local dust storms. 

4.4. Comparison With Other Spacecraft 
Missions 

The Mariner 9 mission arrived at Mars in 1971 dur- 
ing the decay of perhaps the largest dust storm ever 
recorded [Martin and Zurek , 1993]. The Mariner 9 
IRIS instrument [Hanei et a/., 1972a, b] took ther- 
mal infrared spectra similar to TES from L,=293° 
to 352°, with an additional very few spectra between 
L,=352° and 98°. This corresponds in season to TES 
data taken during the science-phasing orbit. Fen- 
ton et al . [1997] convolved IRIS spectra with the 

Viking IRTM spectral response functions and used 
the method of Martin [1986] to retrieve dust opac- 
ities. They found 9 -/mi opacities of 0,5-0. 8 at the 
beginning of the mission soon after the peak of the 
great dust storm (L,=293°-300°). with opacities de- 
caying to about 0.2 after the dust storm had cleared 
(L,=330°-340 o ). Fenton et al. [1997] computed ex- 
ponential decay constants for the dust storm of 42 
days (41 sols) at high southerly latitudes (55°-65°S) 
and 67 days (65 sols) at midsoutherly latitudes (20°- 
30° S). Conrath [1975] found a decay rate of 60 days 
(58 sols) for the 1971 storm based on 2-mbar temper- 
atures. There was no global dust storm observed by 
TES during this season, but the “background” dust 
opacity of 0.2 is similar to that observed by TES (be- 
fore the Z,=309° regional dust storm), and the decay 
constant for the 1971 global dust storm is comparable 
to that for the Noachis storm (49 sols in the core, 75 
sols just to the north) observed by TES. 

Mariner 9 dust opacities have also been retrieved 
by Santee and Crisp [1993], who performed a simulta- 
neous retrieval of atmospheric temperatures and dust 
opacity on IRIS spectra taken between L 5 =343° and 
348° after the dust storm had cleared. They found 9- 
pm dust opacities of 0.25-0.3 near the equator and 
somewhat lower values (0. 1-0.2) to the north and 
south. The variation by a factor of 2 in dust opacity 
with latitude is not observed by TES, but the over- 
all median dust opacity reported by Santee and Crisp 
(0.2) is similar to that observed by TES. 

Dust opacities were monitored by the Viking or- 
biters and landers for more than one Mars year (1976- 


1978). Martin [1986] and Martin and Richardson 
[1993] used data from the infrared thermal mapper 
(IRTM) instrument to retrieve 9-pm opacities, while 
Pollack et al. [1977. 1979] and Colburn et al. [1989] 
used observations of the Sun by the landers to retrieve 
dust opacities at visible wavelengths. TES measures 
infrared opacities, so the IRTM data set is the most 
directly comparable Viking observation. Comparison 
of infrared and visible dust opacities is more difficult 
because the ratio of dust opacity between the two 
wavelengths is not well known. We choose a con- 
version factor of 2.0 for the ratio of visible to 9- pm 
opacities based on the work of Clancy et al. [1995]. 
During the season observed by TES (L J = 184 0 -28 0 ), 
Viking observed two global dust storms, each having 
opacities exceeding unity over a wide latitude range 
[Martin, 1986], with decay constants of 75 sols for the 
first storm and 51 sols for the second storm [Pollack 
et a/., 1979]. The Noachis dust storm observed by 
TES did not reach the size of the storms observed by 
Viking, but the decay constants for all three storms 
were similar. 

Dust opacities observed at Z, s =40° (well away from 
regional dust storms in both the Viking and Global 
Surveyor years) were 0. 1-0.2 using Viking IRTM data 
[Martin, 1986], about 0.5 in the Viking optical data 
[Colburn et al 1989], and 0.1 for TES. Given the 
IRTM spectral response and the dust spectral shape 
from Figure 2, dust opacities from the IRTM should 
be 0.835 times as large as dust opacities reported here 
for TES. Therefore the equivalent IRTM dust opacity 
is 0.12-0.24. Given a measured surface pressure at 
the Viking Lander 1 site of about 8.3 mbar [Hess et 
al. , 1980], and using a conversion factor of 2.0 for the 
ratio of visible to 9-pm opacities [Clancy et al ., 1995], 
the equivalent Viking lander dust opacity becomes 
0.18. Therefore TES dust opacity appears to be some- 
what lower than the Viking opacities, although a new 
calibration of atmospheric temperatures derived by 
the IRTM (R. J. Wilson and M. I. Richardson, The 
Martian atmosphere during the Viking Mission, 1, In- 
frared measurements of atmospheric temperatures re- 
visited, submitted to Icarus , 1999) may affect dust 
opacities derived from IRTM data. At L s — 200°, 
before the major Viking and Global Surveyor dust 
storms, dust opacities were about 0.4 (IRTM dust 
opacity from Mar/ in [1986] converted to TES equiv- 
alent), 0.4 (Viking lander visible dust opacity from 
Colburn et al. [1989] converted to TES equivalent), 
and 0.15 from TES. TES observations imply that dur- 
ing this last Mars year, the southern spring and sum* 
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mer seasons were less dusty than the corresponding 
season during the first Viking year. 

The Mars Pathfinder landed on the surface of Mars 
on July 4. 1997, and operated for 84 sols [Golombek 
et a/., 1997], corresponding to a seasonal range of 
Lj = 142°-189°, during which a short period of over- 
lap exists between TES and Pathfinder dust opacity 
measurements. Smith et al. [1997] used observations 
of the Sun by Pathfinder to obtain dust opacities at 
visible wavelengths. They found dust opacities grad- 
ually rising from about 0.45 at the beginning of the 
mission to 0.55 at the end. To convert to a TES 
equivalent dust opacity, we use a measured surface 
pressure of about 6.7 mbar at the Pathfinder land- 
ing site [Schofield et a/., 1997] and a conversion factor 
of 2.0 for the ratio of visible to 9-^im dust opacities 
[Clancy et a/., 1995]. The TES equivalent Pathfinder 
dust opacities are then 0.20-0.25. This is somewhat 
higher than the 0.15 measured by TES at the latitude 
of the Pathfinder (19°N) and that season. 

5. Summary 

Even before achieving its final mapping orbit, the 
Thermal Emission Spectrometer has already provided 
a very valuable data set by monitoring the Mars at- 
mosphere for more than half a Mars year (L, = 184°- 
28°). During this period, TES observed one large 
regional dust storm (Noachis. at £,=225°), three 
smaller (but still significant) regional dust storms 
(L 5 = 192 0 , 252°, and 309°), and numerous local dust 
storms. During southern spring and early summer 
when the south seasonal polar ice cap was retreat- 
ing, there was consistently more activity (local dust 
storms) and a higher background dust opacity at high 
southerly latitudes near the edge of the ice cap than 
at other latitudes. Nominal dust opacities when there 
were no dust storms were about 0.15-0.2. 

There was no global dust storm of the magnitude 
of those observed by Mariner 9 or Viking observed 
by TES, but the large, regional dust storm (Noachis) 
that was observed by TES was extensive, and it was 
recorded in unprecedented detail. Atmospheric tem- 
peratures and dust opacities were observed before, 
at the height, and during the decay of the Noachis 
storm without interruption. In the sols leading up to 
the Noachis dust storm, there was a significant and 
general rise in southern hemisphere dust opacity and 
in local dust storm activity. The Noachis dust storm 
spread and intensified to its maximum extent in about 
10 sols, covering all longitudes poleward of 40°S lati- 


tude with significant dust opacity (> 0.5). The dust 
storm then decayed with an exponential time constant 
of about 49 sols. A rapid temperature increase of 10- 
15 K at the 0.3-mbar pressure level was observed in 
both hemispheres [Conrath et al this issue] imply- 
ing strong interaction between dust opacity and the 
dynamical state of the atmosphere. 

There is significant and poorly understood interan- 
nual variability in dust opacities on Mars, especially 
in the strength and number of regional and global 
dust storms. The data set returned by TES during 
aerobraking and science phasing helps to extend the 
record of Mars dust observations and also gives us our 
best ever look at the detailed distribution and evolu- 
tion of dust in the Mars atmosphere. 
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Figure 1 . Examples of typical daytime Thermal Emission Spectrometer (TES) spectra. Radiances are shown in 
terms of equivalent brightness temperature to emphasize features. The signature of COo gas, H 2 0 ice aerosols, 
and dust aerosols are all easily identified and are distinct from each other. The uppermost curve shows a hot, 
midday spectrum with nominal dust opacity (0.2). The two lower curves show a spectrum taken at the height of 
the Noachis dust storm and a spectrum taken on the flank of Arsia Mons where there were a thick H 2 0 cloud and 
low dust opacity due to the high elevation. 

Figure 2. The spectral dependence of dust opacity (solid line) normalized to unity at 1075 cm"l. Opacities are 
not computed in the spectral region where C0 2 has significant opacity (560-770 cm ! ). 


Plate 1 . Summary of column-integrated dust opacities observed by TES. Opacities are at 1075 cm 1 and have 
been normalized to an equivalent 6.1 mbar surface. Regional dust storms are apparent at L, = 190°, 225° , and 309°. 

Plate 2. Map showing all retrieved dust opacities taken during the main part of the Noachis dust storm (L,=225°- 
233°). The storm began in the Noachis Terra region around 20°-50°S latitude, 320°-360°W longitude. No quanti- 
tative opacities were obtained at latitudes poleward of 68° S (over the polar cap). The color scale for this and the 
remaining plates represents dust opacity 0.2 as blue, 0.4 as green, and 0.75 or higher as red. Spot size is unrelated 
to the actual footprint size, which ranged from 10 km at periapsis to about 500 km for many of the high southerly 
observations taken in the apoapsis portion of the orbit. 


Plate 3. Orbit-by-orbit evolution of the Noachis dust storm. The orbital period at this time was about 30 hours. 
Large orbit-to-orbit changes in the spatial extent and intensity of the storm were observed. The storm preferentially 
spread to the south and east. The center of activity at 60°~70 o S and 180° -240° W persisted for quite some time 
after the main center in Noachis had largely subsided. 

Plate 4. Map showing all retrieved dust opacities taken during the regional dust storm that occurred during 
science-phasing orbit (L,=309°-3 19°). The storm was nearly as intense as the Noachis storm but was much more 
confined spatially. 

Plate 5. Map showing all retrieved dust opacities taken during the regional dust storm that occurred at the 
beginning of the mission (L, = 184°-195°). Two centers of activity are seen near Hellas and Argyre, along with a 
general increase in dust opacity along the edge of the retreating south polar seasonal ice cap. 
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Table 1 . Significant Local Dust Storms Observed by TES 


L,, deg 

Location 

Comments 

205 

25°S, 295° W 

small; observed in only one orbit 

208 

40-60°S, 170-200®W 

small; observed in two orbits 

219 

40-60°S, 1 00-1 80° W 

beginning of numerous widespread small 
storms along cap edge 

224 

45°S, 130°W 

simultaneous with birth of Noachis storm; 
observed in only one spectrum 

236 

30°N, 135° W 

high opacity for north hemisphere 

252 

35-45° S, 1 30-170° VV 

large, intense storm but gone quickly 

257 

70°S, 350-360® W 

small; beginning of much activity at high 
southerly latitudes 

258 

70-75° S, 165-210°W 

related to Noachis satellite? intense and 
lasted several orbits 

261 

70-80°S, 250-290° W 

spread east and intensified 

267 

70-80°S, 190-230® W 

intensification there part of a general 
dustiness at high southerly latitudes 

304 

20® N, 45® VV 

near Viking 1 lander; also observed with 
the MGS MOC camera 

0 

20-30° S, 330-360®W 

late in the season; died away quickly 

3 

40-50°S, 40-60® W 

late in the season; died away quickly 
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